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Abstract 



Interacting white dwarf binary star systems, including helium cataclysmic 



variable (HeCV) systems, are expected to be strong sources of gravitational 

psj ■ radiation, and should be detectable by proposed space-based laser interfer- 

on ; 

ometer gravitational wave observatories such as LISA. Several HeCV star 
systems are presently known and can be studied optically, which will allow 
5— i . electromagnetic and gravitational wave observations to be correlated. Com- 

parisons of the phases of a gravitational wave signal and the orbital light curve 
from an interacting binary white dwarf star system can be used to bound the 
mass of the graviton. Observations of typical HeCV systems by LISA could 
potentially yield an upper bound on the inverse mass of the graviton as strong 
as h/nig = X g > 1 x 10 15 km (m g < 1 x 1CP 24 eV), more than two orders of 
magnitude better than present solar system derived bounds. 
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I. INTRODUCTION 



The advent of large scale laser interferometer gravitational wave detectors promises to 
lead to the direct detection of gravitational radiation, which will create a new observational 
field of science: gravitational wave astronomy. The introduction of gravitational waves 
into the retinue of astronomical observations will yield important new information about 
the dynamics of astrophysical systems, and will also provide an excellent opportunity to 
conduct new tests of gravity. 

Proposed space-based laser interferometers, such as the LISA (Laser Interferometer Space 
Antenna) [|IJ and OMEGA (Orbiting Medium Explorer for Gravitational Astrophysics) || 
observatories, will be particularly well poised to begin astrophysical studies since there are 
known sources of gravitational radiation which will be easily visible to these instruments, 
namely interacting binary white dwarf (IBWD) star systems. 

The IBWD sources are particularly appealing targets because they are good candidates 
for simultaneous optical and gravitational wave observations. Many of the IBWDs which are 
known to exist are being studied and monitored by the Center for Backyard Astrophysics 
(CBA)[|, and are expected to be strong sources of monochromatic gravitational waves which 
should be easily visible to an instrument such as LISA with only a few minutes of sig- 
nal integration. Simultaneous optical and gravitational wave observations will be useful in 
refining the current physical models used to describe these systems, and for testing rela- 
tivistic theories of gravity in the radiative regime by comparing the propagation speeds of 
electromagnetic and gravitational wave signals. 

This paper examines how the comparison of the phase of the orbitally modulated electro- 
magnetic signal (the light curve) and a gravitational wave signal from an IBWD star system 



The CBA is a network of amateur astronomers equipped with CCD photometry equipment 
who monitor variable stars. The network is managed by professional astronomers at Columbia 
University. 
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can be used to bound the mass of the graviton. If the mass of the graviton is assumed 
to be known by other measurements, then the observations may be used to determine the 
properties of the binary star system being monitored. 

Current conservative bounds on the graviton mass come from looking for violations of 
Newtonian gravity in surveys of planetary motions in the solar system. If gravity were 
described by a massive field, the Newtonian potential would have Yukawa modifications of 
the form 

V(r) = ~exp(-r/X g ) (1) 

where M is the mass of the source of the potential, and X g = h/m g is the Compton wave- 
length of the graviton, where m g is the graviton mass. The current best bound on the 
graviton mass from planetary motion surveys is obtained by using Kepler's third law to 
compare the orbits of Earth and Mars, yielding X g > 2.8 x 10 12 km (m g < 4.4 x 10~ 22 eV) 

!• 

Another bound on the graviton mass can be established by considering the motions of 
galaxies in bound clusters Q], yielding \ g > 6 x 10 19 km (m g < 2 x 10~ 29 eV). This bound, 
while stronger than solar system estimates, is considerably less robust, due to uncertainty 
about the matter content of the Universe on large scales (e.g., the amount and nature of 
dark matter is widely debated, and uncertain at best). 

Recent work by Will f5[ has suggested that the mass of the graviton could be bounded 
using gravitational wave observations. If the graviton is a massive particle, then the speed of 
propagation of a gravitational wave will depend on its frequency. As binary systems evolve, 
they will slowly spiral together due to the emission of gravitational radiation. Over the 
course of time, the frequency of the binary orbit rises, ramping up rapidly in the late stages 
of the evolution, just prior to coalescence. Laser interferometer gravitational wave detectors 
should be able to track the binary system's evolution, obtaining the detailed time-dependent 
waveform using the matched filtering techniques required for data analysis in these detectors. 

Space-based detectors such as LISA will be able to observe the coalescence of massive (~ 
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10 5 to ~ 1O 7 M ) binary black holes, as well as the gravitational wave emission from compact 
binary star systems which are far from coalescence (e.g., interacting binary white dwarfs). 
Ground-based detectors such as LIGO will be able to detect the merger of smaller black 
hole binaries (~ 10M Q ), as well as the coalescence of compact binary stars (e.g., neutron 
star/neutron star binaries). If the graviton is a massive particle, then the observed signal will 
not perfectly match theoretical templates computed using general relativity theory, in which 
the graviton is massless; a massive graviton would cause dispersion in the gravitational waves. 
By using matched filtering of inspiral waveforms, this dispersion could be bounded, thereby 
bounding the mass of the graviton. Will || finds that LIGO could bound the graviton mass 
at X g > 6.0 x 10 12 km (m g < 2.1 x 1CT 22 eV) by observation of the inspiral of two 10 M 
black holes. A space-based interferometer such as LISA, observing the inspiral of two 10 7 M Q 
black holes could bound the graviton mass at X g > 6.9 x 10 16 km (m g < 1.8 x 10~ 26 eV). 
If the graviton is massive, then these numbers represent the minimum masses detectable by 
such observations. 

The analysis in this paper shows that LISA observations of known IBWD sources could 
yield a bound as strong as X g > 5 x 10 15 km (m g < 2 x 10~ 25 eV), considerably stronger 
than present solar system based bounds. The IBWD bound also has the advantage of not 
depending on the complicated details of black hole coalescences. 

Section reviews what is known about the interacting binary white dwarfs, in particular 
the helium cataclysmic variable (HeCV) systems and their archetype, the binary AM CVn 
(AM Canum Venaticorurri) . Section |T| reviews the basic notions associated with (possibly) 
massive photons and gravitons. Sections [TV]— [V] propose a new experiment to measure the 
graviton mass using IBWD observations, and an expression for the mass is derived. In 
Section |VT|, the sensitivity predicted for LISA is used to estimate how precise a bound could 
be placed on the graviton mass from HeCV observations. Section |V11| summarizes the results, 



and also suggests how correlation of phase measurements might be used to measure other 
astrophysical parameters in the binary system, such as the accretion disk radius. 

Throughout this paper, geometric units with G = c = 1 are employed, unless otherwise 



noted. 



II. INTERACTING BINARY WHITE DWARF SYSTEMS 

Estimates suggest the Galaxy is populated by ~ 10 7 close binary star systems ||. The 
sheer numbers of these systems is likely to have profound consequences for space-based 
gravitational wave observatories. The combined gravitational waves from these binaries 
will produce a stochastic background which rises well above the low frequency detection 
sensitivity of LISA |7| . Particularly strong (e.g., nearby) binary systems will rise above 
this background and be observable by a spaceborne observatory. One class of such sources 
are the helium cataclysmic variable (HeCV) stars. The properties of the six nearest known 
HeCVs are shown in Table p] ||. The predicted stochastic gravitational wave background 
due to short period binary stars, as calculated by Hils and Bender J?| and the predicted 
signals of the six nearest HeCVs are plotted in Figure [I], along with the predicted sensitivity 
curve of LISA. If present models of the spatial density of close binaries in the Galaxy are 
correct, roughly 5000 of these sources should be individually detectable by a space-based 
laser interferometer such as LISA Q. 

Currently the best models for HeCVs describe a star system where the secondary star 
(the lower mass companion in the binary system, usually a degenerate helium dwarf star) has 
expanded to fill its Roche lobe, and the primary star (the larger mass, compact white dwarf) 
lies at the core of an accretion disk. Matter overflows from the secondary Roche lobe and 
streams onto the accretion disk, creating a hot spot which emits a strong electromagnetic 
signal. 

There are several mechanisms whereby the light curve of an HeCV could be modulated 
as seen from the perspective of observers on Earth. The simplest model is for systems whose 
orbital plane is close to the line of sight to the Earth, so that the stars periodically eclipse, 
partially or wholly. The eclipse phase could dim the stellar components, part of the primary 
emission regions of the accretion disk, or the Roche lobe. Another possible mechanism for 
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variation in the light curve is associated with the emission from the secondary star, which has 
expanded to fill its Roche lobe. The star will appear brighter when the largest surface area 
is presented to the observer along the projected line of sight. This will occur twice in each 
orbit, when the line of sight is perpendicular to the line of centers. Signatures in the light 
curve due to projected area effects are called "ellipsoidal variations," and are most easily 
observed in infrared wavelengths |J. Yet another possible mechanism is currently favored 
to explain the source of variation in the light curve of the archetype system for HeCVs, 
AM CVn [ 10 1 . In this model the hot spot on the accretion disk radiates approximately 



radially outward from the disk. As the binary orbits, this hot spot alternately turns towards 
and away from distant observers, leading to a modulation of the light curve (a so-called 
"flashlight" mechanism). 

A detailed theoretical model of AM CVn has been constructed, describing a variety of 
signals which are present in the photometric data JTTHT3|. This model suggests that AM 
CVn is a member of a class of variable stars that have periodic features in the light curve 
known as "superhumps" [ [TjjTS ]. The model explains the superhump feature as being caused 



by the existence of an eccentric precessing accretion disk, with a precession period which is 
slightly longer than the orbital period of the binary. Knowing the superhump period, P s h, 
and the precessional period of the accretion disk (apsidal advance), P aa , the model predicts 
the orbital period will be given by 

Porb = P s h + P<ia ■ (2) 

Photometry of AM CVn shows the existence of a superhump signature at P s h = 1051.2 s, 
and the period of the accretion disk precession at P aa = 13.38 hr. Using Eq. (|2|) this model 
predicts a binary orbital period of P or f, = 1028.77 ±0.18 s for AM CVn |Tj|. Photometric 
observations by the CBA have recently confirmed an orbital period of P or b = 1028.7325 ± 
0.0004 s 0. 
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III. MASSIVE GRAVITONS AND PHOTONS 



To constrain the mass of the graviton by comparing the propagation speed of gravita- 
tional and electromagnetic waves, one must consider how the speed of gravitational waves 
(and electromagnetic waves) is related to the mass of the graviton (and the possible mass 
of the photon) . 

The current bound on the photon mass is m 7 < 2 x 10~ 16 eV |[L6|| , which is much larger 
than the current bounds on the graviton mass (cf., the solar system bound on the graviton 
mass is m g < 4.4 x 10 -22 eV). Is it then justifiable to treat the photon as a massless particle, 
while at the same time treating the graviton as a massive particle? 

The resolution to this question can be understood by examining the partition of energy 
between the rest mass and kinetic energy of a particle being received from a distant binary. 
From the relativistic energy, E 2 = p 2 + m 2 , one may write the velocity of any particle as 

(3) 



If m <C E, then Eq. (|) implies 



1 - v ~ , (4) 

IE 2 K J 



where e parameterizes the difference between the velocity of the particle and c. 

For optical photons (A ~ 500 nm) received from a binary star system, the characteristic 
energy is E 1 ~ 2.5 eV. For this energy, e 7 < 3 x 10~ 33 . Similar considerations may be 
applied to the gravitons received from the same binary systems. In this case, the frequency 
of the gravitational waves is / ~ 10~ 3 Hz, giving a characteristic energy for a single graviton 
of Eg ~ 4 x 10~ 18 eV. Using the solar system bound on the graviton mass, m g < 4.4 x 10~ 22 
eV, yields e g < 1 x 10~ 8 . 

For the current bounds on the photon and graviton masses, e 7 <C e g . If the bound on 
the mass of the graviton is not drastically improved (e.g., decreasing the bound on m g by 
12 orders of magnitude, such that e g ~ e 7 ), then the effect of a non-zero mass will be much 
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more significant for gravitons than photons in our analysis. This justifies the treatment of 
the photons as massless particles and the gravitons as massive particles in this paper. 

Since the gravitational waves emitted by an HeCV binary star system are essentially 
monochromatic, they will have a single particular velocity, 



Writing the energy in Eq. (j5|) as E = hu, and identifying the Compton wavelength of the 
graviton as X g = h/m g this becomes 



IV. CORRELATION OF ELECTROMAGNETIC AND GRAVITATIONAL 
OBSERVATIONS TO MEASURE THE GRAVITON MASS 

As noted in Section |T|, it is expected that the Galaxy harbors a large population of 
interacting binary white dwarf stars. The light curves for several of these systems are 
already known, obtained from ground-based optical photometry. Because these systems are 
expected to be observable in the gravitational wave spectrum (in addition to the optical 
spectrum) , they present an excellent opportunity to directly compare the propagation speed 
of electromagnetic and gravitational waves. 

Consider the schematic diagram shown in Figure 0. The phase fronts of the light curve 
modulation are represented in the top half of the diagram. The binary star system will also 
emit gravitational radiation which could be monitored by Earth-bound observers as well. 
The phase fronts of the gravitational wave signal are represented in the lower half of the 



2 Assuming a circular orbit, the frequency of the gravitational radiation is actually twice the orbital 
frequency of the binary. For clarity, only half the gravitational phase fronts have been drawn in 
the diagram. 
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diagram^]. 
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Suppose the two signals are emitted in phase at the source as shown. If the graviton 
is a massive particle, then the gravitational waves propagate at a speed v g < 1, and the 
gravitational phase fronts will lag behind the light curve phase fronts when the signals 
arrive at Earth, as shown. By measuring the lag between the phase fronts, the mass of the 
graviton can be measured or bounded. 

To determine the lag between the two signals, the phase of each signal must be measured. 
Consider a binary with orbital frequency u at a distance D from Earth. Assuming the 
photon to be massless (v = c — 1), the observed phase of the light curve will be 

(j) em = uj D + A , (7) 

where the term A represents a variety of effects (discussed in Section [V|) which could create 
phase delays between the electromagnetic and gravitational signals that are being monitored. 
In contrast, the gravitational wave signal, traveling at v g < 1 will arrive at Earth with a 
phase of 

4>gw = 2^ — . (8) 
V 9 

The phase lag, A, between the light curve and gravitational wave signals is constructed 
from these two phases: 

= oj D ( 1 - 1 + a) , (9) 

where a = A/(u D). The factor of 1/2 insures that the phase subtraction is done between 
two signals with the same frequency. It is convenient to define the fractional change in the 
phase as 

e=-^-= (1-1 + a) . (10) 

Taking the definition of v g from Eq. ([]) and substituting into Eq. (|T0"D, the Compton 
wavelength of the graviton as a function of the fractional phase lag is found to be 
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— a 



1 



) 



(11) 



An obvious question to ask about this analysis is how to determine whether or not the 
phase difference between the two signals is greater than a single cycle, and hence unde- 
tectably large (e.g., A = 2mr + £, where £ is a small quantity). One can eliminate this 
concern because strong bounds on the graviton mass already exist. The largest observable 
phase shift which is consistent with current bounds can be computed by simply evaluating 
Eq. ([11]) with X g equal to the bound of interest. For example, the bound on the graviton 
mass given by solar system constraints, applied to the AM CVn system, yields a maximum 
fractional phase change of e = 1.5 x 10 -9 . Since AM CVn lies at a distance of D — 101 pc, 
this value of e indicates a maximum phase difference of A = 9.6 x 10~ 2 for that system (we 
have let a — > here for convenience; if the measured phase difference were larger than this 
value, it would indicate that q^O). 



In order to evaluate Eq. (p]), one must not only measure the phase lag between the two 
signals, but an estimate must be made for the value of a. 

The parameter a can be written as the sum of two primary sources of delay between the 
gravitational and electromagnetic signal phases: 



where a pat h is a phase lag associated with the wave's propagation from the binary to the 
observer at the Earth and a* is a phase lag which depends on the specific astrophysical 
nature of the binary star system. 

In principle, a pat h will be nonzero because the line of sight to the binary is an imperfect 
vacuum, with non-unit index of refraction. The variations in index of refraction over the 
path will cause a lag in the electromagnetic signal. The dominant source of this lag will be 
caused by propagation of the signal through the Earth's atmosphere. 



V. PHASE DELAYS 



a = a* + ap a th , 



(12) 
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A simple estimate of the value of a pa th can be made by computing the electromagnetic 
phase delay due to propagation through a modeled exponential atmosphere, with a density 
profile 

p(r) = p exp(-r//i s ) , (13) 

where h s is the scale height of the atmosphere. If the index of refraction, n, is assumed to 
vary linearly with the density, then 

n(r) = 1 + ?7exp(— r/h s ) , (14) 

where i] = (n atm — n vac ) = (n atm — 1) is the difference in index of refraction between the 
atmosphere and vacuum. The index of refraction is related to the signal's propagation speed 
v by 

dr 1 



dt n(r) 

Eq. ( |i~5D can be integrated using Eq. (|T^) to obtain 



(15) 



\ ttmi dt= r° dr[l + r]exp(-r/h s )] , (16) 
Jo Jo 

where t tra nsit is the time it takes a photon to transit the atmosphere, and r Q is the height at 
which the effects of the atmosphere become negligible. Completing the integration yields 

ttransit = r Q + h s 7]{l - e~ r ° /hs ) . (17) 

In order to compute a phase delay, one is interested in the time by which the photons are 
delayed by the atmosphere, which is 

tdelay ^transit ^vacuum 

^transit '"o > (1^) 

where t vacuum is the time it would take a photon to travel the same distance (i.e., the 
atmospheric depth) in vacuum. 
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The phase delay introduced by this effect is obtained by multiplying the delay time by 
the frequency of the signal being observed, 

4> delay = ^ot delay ■ (19) 

The parameter a pat h is obtained from the phase delay by dividing by the reference phase, 
lo D (as in Eq. (g)), giving 



— ^ = ^- < 20 > 

It is now possible to numerically estimate the delay introduced by the atmosphere. 
Assuming the Earth's atmosphere to be in hydrostatic equilibrium gives a scale height 
h s = 8500 m. A typical value for the difference between the atmosphere's index of re- 
fraction and unity, evaluated at sea level, is 2.8 x 10 -4 . Taking the atmospheric depth to 
be r Q = 10h s , Eq. ( |T8| ) gives t de i ay = 8.0 x 10~ 9 s. For the prototypical source AM CVn, 
at a distance D = 101 pc, Eq. fl20"|) yields a pat h = 7.7 x 10 -19 . This is much smaller than 
any possible phase measurement, either electromagnetic or gravitational, and so we will 
henceforth ignore this source of uncertainty. 

The parameter a* is a measure of the initial phase difference at the source between the 
gravitational wave and electromagnetic signals. It indicates the relative phase difference 
between the peaks in the light curve of the binary, and the peaks in the quadrupole gravi- 
tational radiation pattern. Determining the value of a* requires knowledge of the position 
of the stars in the binary system when the electromagnetic signal peaks. The quadrupole 
gravitational radiation pattern will peak along the line of masses in the binary system, and 
also 180° away from the line of masses (since the frequency of the gravitational radiation is 
u = 2uj ). If the primary variation in the binary light curve is associated with the transit of 
the line of masses across the observer's line of sight (e.g., the system is an eclipsing binary), 
then it is straightforward to assign a+ = 0, indicating no initial phase delay between the 
gravitational wave signal and the binary's light curve. 

For more complicated systems such as AM CVn, the light curve variation reflects the 
orbital motion of the hot spot on the edge of the accretion disk where the matter stream 
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from the secondary Roche lobe overflow strikes the disk. Studies of Roche lobe overflow 
fT7|| show that the transferred material remains in a coherent matter stream which spirals in 
towards the primary star. In this sort of system, the location of the hot spot will determine 
the value of a*, which will describe the amount by which the hot spot leads the line of 
masses, as shown in Figure |3|. 

Estimates of the size of the primary accretion disk in HeCV type systems suggest that 
disk radii will be around 75% of the primary Roche radius |18|], but this estimate is only 



certain to within about 10%. This uncertainty makes it virtually impossible to estimate a* 
adequately for use in Eq. (|TT|) from present observational data . Future observations of HeCV 
systems, either from advanced ground-based instruments such as the Keck Interferometer, 
or from space-based instruments such as the Space Interferometry Mission ]19 and the 



Terrestrial Planet Finder |2(| could allow direct measurement of a* by optically imaging the 



detailed structure of close binary systems. 

For cases when a* cannot be accurately determined, the dependence of the Compton 
wavelength of the graviton on this parameter (or on a pat h) can be eliminated by subtraction 
of two observations of the source. Consider the situation shown in Figure [|, where the 
gravitational and electromagnetic signals are monitored when the Earth lies on one side of 
its orbit, and again six months later, when it lies on the other side of its orbit. 

When the Earth is in position 1, the phase difference between the electromagnetic and 
gravitational wave signal can be written 

uJoD _ UoD + A (21) 

Similarly, when the Earth is in position 2, the phase difference may be written 

a u>o(D + L) 
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u (D + L) + A, (22) 



where L < 2AU is the path length difference for the time of flight between the two measure- 
ment. Subtraction eliminates the unknown quantity A, yielding 

A 2 - A 1 = u l(^--1\ . (23) 
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Defining the fractional change in phase from this quantity gives 



e = 



A 2 - A 




(24) 



which in terms of the Compton wavelength becomes 




(25) 



The unknown parameter, a, has been eliminated from the expression for A 9 , but at the cost 
using a much shorter characteristic distance, L <^ D. This approach amounts to measuring 
the phase lag between the two signals over the time it takes to cross the Earth's orbit (at 



In order to estimate a bound on the graviton mass, assume a null result for the measure- 
ment of the phase difference, A, between the two signals. The size of A (and therefore e) 
will then be limited only by the uncertainty in the measurements of the phase. Combining 
the uncertainty of the gravitational phase measurements with the electromagnetic phase 
measurements in quadrature yields 



where S(f) em and S(j) gw are the uncertainties in each of the phases. 

For observations with a space-based interferometer, the error in phase measurements can 
be estimated as the ratio between the sampling time and the total integration time. For 

3 An improvement of roughly a factor of two may be obtained by combining the Earth's orbital 
motion with the proper motion of the binary system relative to the Sun, which will typically be of 
order a few AU per year 



most)[], as opposed to the time it takes to propagate over the Earth-source distance, leading 



to a great loss of precision in the measurement of X g . 



VI. OBTAINING A MASS BOUND 




(26) 
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LISA, the sampling time is expected to be of order 1 s, with total integration times of 1 yr 
= 3 x 10 7 s, yielding 8(f> gw = 3 x 10 -8 . The CBA reports a 0.0004s uncertainty over the 
1028.7325s period of AM CVn, yielding a phase uncertainty of 5<j) em = 4 x 10~ 7 . 

For the case of AM CVn, the value of a* is still not known, so bounds on the graviton 
mass must be derived from Eq. fl2"5]). As shown in Figure [|, the characteristic distance L 
is simply the path length difference between the two measurements. If the inclination of 
the binary system to the plane of the Earth's orbit is j3, then the characteristic distance 
is L = 2cos(/9)AU. For AM CVn, which lies at ecliptic latitude (3 = 37.4°, this yields 
L = 2.38 x 10 11 m. Using this value, Eq. (|2~5"D gives a bound 

A s >5x 10 14 m = 5 x 10 11 km , (27) 

or, in terms of the graviton mass, m g < 2 x 10~ 21 eV, about a factor of five worse than the 
present bound based on the motion of Mars. Even this weak bound would be of interest, 
however, since it is based on the dynamics of the gravitational field, (i.e., gravitational 
waves, rather than the static Yukawa modifications of the Newtonian potential). 

If the value of a+ could be determined precisely (e.g., by monitoring ellipsoidal variations 
the light curve in the infrared, as suggested in Section [TT], or with future optical interferometer 
observations), such that the uncertainties 5a± ^ 10~ 7 , then Eq. could be used to bound 
the graviton mass. The distance to the known IBWD systems is typically of order D ~ 100 
pc; combining this with a typical orbital period of P ~ 1500s yields a bound 

A s >lx 10 18 m = 1 x 10 15 km , (28) 

or m g < 1 x 10~ 24 eV. This potential bound would be a factor of four hundred more stringent 
that the present solar system based bound, and would be better than the bounds obtained 
from inspiraling black holes proposed by Will || for all but very large black holes. 

During the next decade, as we await the launch of LISA, the optical astronomers may 
well succeed in further reducing the uncertainty in their phase measurements. If the optical 
signal phase error is reduced in Eq. (]2q) to the point where the dominant source of error is the 
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gravitational wave phase measurement, then the bound obtained above could be improved 
by about another factor of five, to X g > 5 x 10 15 km (m g < 2 x 10~ 25 eV). 



VII. SUMMARY 

After the initial detection of gravitational waves the challenge will be for the field to 
evolve into a productive observational science which makes firm contact with astrophysics, 
complementing the broad base of electromagnetic observations already supporting that field. 
The experiment proposed here is particularly appealing because it entails observations of 
known sources by space-based detectors. The existence of IBWDs has been verified (as 
opposed to more speculative sources, such as binary black hole coalescence events), and 
such objects are currently under study by observational astronomers. Detailed gravitational 
wave observations can begin almost as soon as a space-based interferometer such as LISA is 
online. 

We have shown that reliable bounds on the mass of the graviton of order X g > 1 x 10 15 
km could be obtained through detailed observations of the interacting binary white dwarf 
star systems such as AM CVn. With the combination of detailed studies of such binary 
systems by optical interferometers and gravitational wave observations, this could be a very 
robust bound, several orders of magnitude greater than the current best bounds from solar 
system observations. 

If one assumes the graviton to be a massless particle, as predicted by general relativity, 
then the same measurements described here can be employed to determine the structure 
of the binary star system. If the graviton is massless, then any phase difference measured 
between the gravitational wave and electromagnetic signal must be due to effects in the 
binary system. Setting v g — 1 in Eq. ( |T0l) yields 

e = ^ = a - < 29) 

showing that the difference in phase is simply an indicator of the value of a = a pa th + As 
was shown in Section |\T|, the value of a pa th is expected to be negligible (a pat h = 7.9 x 1CT 19 for 
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AM CVn). In the cases where a pat h can be ignored, the measured phase difference will be a 
direct measure of the value of A, which is the amount of phase by which the electromagnetic 
signal leads the line of masses in the binary system. With good models of the matter stream 
from the secondary Roche lobe overflow (such as the trajectories shown in Figure |3|), a 
measurement such as this could allow an accurate determination of the accretion disk radius 
and the refinement of physical models for HeCV type stars solely from gravitational wave 
observations. 
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TABLES 

TABLE I. Properties of the six nearest interacting binary white dwarfs. 



Name 


Secondary Mass 
(M Q ) 


Orbital Period 

(s) 


Distance 

(PC) 


Strain h 
(xlO^ 22 ) 


AM CVn 


0.044 


1028.73 


101 


5.27 


EC15330-1403 


~ 0.04 


1119 


165 


5.36 


CR Boo 


~ 0.03 


1491 


157 


2.82 


V803 Cen 


~ 0.03 


1611 


~ 100 


4.20 


CP Eri 


~ 0.03 


1724 


~ 100 


4.02 


GP Com 


~ 0.02 


2791.2 


165 


1.77 
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FIGURES 

FIG. 1. The background of close binaries in the Galaxy, plotted with the projected sensitivity of 
LISA. The six most well understood AM CVn type binaries are indicated. The assumed bandwidth 
is 1(T 7 Hz. 

FIG. 2. Schematic of a binary star system which has observable electromagnetic and gravita- 
tional wave signals. If the graviton is a massive particle, then the phase fronts of the gravitational 
signal will lag behind those of the electromagnetic signal. 

FIG. 3. The Roche lobe of the primary star is shown as a dark line; matter overflows from 
the secondary star in a coherent matter stream, shown as trajectories falling towards the white 
dwarf (indicated by * at the center of the Roche lobe; figure adapted from Flannery(1975) ). 
Depending on the radius of the accretion disk, the matter stream will strike the disk at some angle 
9 which leads the binary axis; this angle is directly proportional to the phase delay parameter a*. 
The dashed circle represents an accretion disk radius which is 70% of the Roche radius. A 10% 
uncertainty in the radius of the accretin disk will yield ~ 5° uncertainty in the value of 6. 

FIG. 4. Schematic showing two observations of a source from opposite sides of the Earth's 
orbit. The radius of the Earth's orbit is te and the path length difference for signals propagating 
from a distant source is L. 
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